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1 .  Introduction 


The  aim  of  this  project  is  to  investigate  vortex-profile  interaction.  Vor¬ 
tex  interaction  with  solid  walls  generates  sound.  Vortex  noise  occurs  near¬ 
ly  everywhere  in  aerodynamics  and  especially  it  is  one  of  the  helicopter 
noise  sources. 

There  is  lack  of  knowledge  about  the  near-field  flow  at  the  profile. 
Theoretical  considerations  are  difficult,  because  the  role  of  the  Kutta- 
Condition  is  not  understood  and  the  full  Navier-Stokes  equations  are  not 
easy  to  solve.  Thus  experiments  must  be  carried  out.  After  some  consider¬ 
ations  about  vortex  paths  with  the  method  of  conformal  mapping,  the  first 
step  into  the  problem  was  an  experiment  which  gave  an  overview  on  the 
phenomena  occurring.  In  a  transonic  tunnel  a  profile  was  placed  in  the  wake 
of  a  vortex  shedding  cylinder.  The  interaction  of  the  Karman  vortex  street 
and  the  profile  can  be  measured  by  pressure  transducers  and  high  speed 
interferograms.  Pressure  waves  and  flow  pattern  can  be  evaluated  from 
these  experiments. 

Additionally,  for  single  vortex  interaction  a  shock  tube  experiment  was 
performed.  Single  starting  vortices  are  generated  with  an  additional  air¬ 
foil.  The  boundary  layer  is  very  thin  in  this  case.  So  the  vortices  are 
similar  to  potential  vortices.  These  experiments  allow  comparisons  to  the 
theoretical  results  achieved  with  the  method  of  conformal  mapping. 
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by  vortex-profile  interaction 


Following  the  work  of  Obermeier,  the  path  of  a  cylinder  vortex  around  a 
Joukowsky  profile  is  calculated  according  to  potential  theory  without  tak¬ 
ing  into  account  the  Kutta  condition.  Sound  generation  is  calculated  by 
matched  asymptotic  expansions  [1,2,3). 


The  flow  potential  of  vortex-profile  interaction 


Potential  flow  of  a  vortex  near  a  wall  is  described  by  two  mirror  vortices 
in  the  complex  u-plane.  The  plane  of  symmetry  corresponds  to  the  wall. 
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Using  the  rational  function 


u+1 

v+1 

— 

u=  - 

u-1 

v-T 

the  imaginary  axis  is  mapped  to  the  unit  circle  in  the  v-plane.  One  gets 
the  potential  of  flow  of  a  vortex  near  a  cylinder 


V 


A  free  stream  UQ  and  an  additional  value  of  circulation  around  the  origin 
are  added  to  this  potential;  circulation  around  the  cylinder  is  thus  eli¬ 
minated  and  the  Kutta  condition  is  satisfied  when  the  vortex  is  far  away. 
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As  a  preparation  for  the  Joukowsky  map,  the  radius  of  the  circle  is  en¬ 
larged  by  a  factor  (a+b)  and  the  center  is  translated  by  an  amount  a.  Tbe 
width  of  the  profile  depends  on  the  parameter  a,  the  roundness  of  the 
trailing  edge  on  b. 
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Using  the  Joukowsky  transform 


Z=  4(w  +  u) 
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one  gets  a  wing  profile  which  approximates  standard  symmetrical  profiles 
except  for  its  rounded  trailing  edge  when  the  parameters  a  and  b  are  chosen 
appropriately.  For  the  complex  root  the  convention 


Re(V  1  -  —  0 

zc 


is  used  m  the  following. 


New  functions  A,  B  and  C  are  introduced  to  simplify  notation: 


B(z )+2a+b-1 


B(z)-b-1 


B(z)+a-1 


-dr  A(z)  -  A(z  )  1 

-  (  i0g  - 2_  +  i0g  C(z)  )  +  H  (C(z)  +  -  ) 

2  n  A(z)  +  A(T0)  °  C(z) 


Note  that  the  velocity  at  infinity  now  is  given  by  L'w  =  2  U0  /(a+b). 

As  a  next  step  the-  nond intensions]  v.u  lable  z,  7./1  j  is  introduced.  (The 
index  l  will  be  su|.pr v.v.r.i<l  for  th<-  :  «Ke  of  simplicity )  .  .The  length  of  the 
profile  equals 

1  =  a+b  +  ( a+b ) / (  ( b+ 1  ) l 2a+b- 1  )  ) . 

Introducing  the  parameter 
c  : =  1  +  1 /  (  (b+1 ) ( 2  a  +b— 1 )  ) 
und  using  the  coefficient  of  circulation 
Z  :=  2  T /(U*  1) 

we  get  the  nondimensional  potential 
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Tfie  potential  of  a  vortex  near  a  profile  may  be  written  in  a  somewhat  dif¬ 
ferent  form: 
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which  is  achieved  by  mapping  the  cylinder  potential 
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Besides  this  complex  representation,  a  polar  coordinate  representation 
using  real  functions  has  been  described  in  the  literature.  Walton  [4] 
used  such  a  representation  for  the  case  of  a  cylinder  and  Lin  [  5  ]  for  a 
flat,  plate. 


The  potential  used  in  this  work  does  not  satisfy  the  Kutta  condition  when 
the  vortex  is  near  the  profile.  This  deficiency  can  be  corrected  by  match¬ 
ing  the  inner  circulation  and  allowing  vortices  to  be  generated  at  the 
trailing  edge.  However,  this  method  encounters  difficulties  because  chaotic 
behaviour  may  appear  when  many  free  vortices  are  present. 


Parthasarathy  [6]  avoids  this  difficulty  by  imposing  restricting  condi¬ 
tions:  the  vortices  generated  at  the  trailing  edge  swim  off  horizontally 
and  with  the  velocity  of  the  surrounding  fluid. 


As  the  author  does  not  know  of  a  method  not  imposing  restricting  condi¬ 
tions,  the  Kutta  condition  has  not  been  taken  into  account  in  this  work. 
Ibis  approach  may  be  justified  as  it  is  not  known  whether,  and  if  so,  in 
which  way  the  Kutta  condition  is  satisfied  in  a  high  frequency  regime,  for 
example  when  vortices  are  passing  nearby. 
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The  path  of  a  vortex 

The  movement  of  a  vortex  can  be  deduced  from  the  translational  velocity  at 
the  location  of  the  vortex  core.  Thus  velocity  is  calculated  by  a  limit 
process:  the  rotational  velocity  near  the  vortex  core  is  subtracted  from 
the  total  velocity.  The  following  expressions  are  valid  for  the  conjugated 
velocity  of  the  core: 
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The  derivatives  of  the  functions  A,  B  and  C  are  given  by: 


A(z)  =  - -  c(z)  = 

(B(2)-b-if,.\ATT 
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An  elegant  method  of  calculation  uses  a  Hamiltonian 
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The  function  H  is  a  constant  of  motion  [3,5,7].  Ibis  facilitates  the  con¬ 
trol  of  accuracy  when  evaluating  the  differential  equation  numerically. 
Evaluation  is  done  by  the  Runge-Kutta  method.  Figure  la  shows  6  vortex 
paths  around  a  flat  plate.  The  coefficient  of  circulation  is  10,  the  vor¬ 
tices  are  rotating  counterclockwise.  The  length  of  the  strokes  corresponds 
to  the  unit  of  time.  There  are  closed  paths  around  the  profile  and  beneath 
it.  Other  calculations  [8,9]  show  that  vortex  trapping  is  possible,  too. 


Tbe  far  field 


When  z  =  z  (t)  is  known,  the  hydrodynamic  far  field  can  be  calculated: 
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Let  T  :=  —5 


t  be  the  nondimensional  time  and 


if 


p.  ;  *  p  -  the  nondimensional  pressure. 

1  rp 


obeys  the  relation  |V  <j>|  ^  ) 


Let  W  be  the  stream  function  of  the  far  field  approximation  of 
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TVien  it  follows  for  the  derivative 
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Furthermore , 
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is  valid. 


The  pressure  in  the  far  field  can  be  expressed  in  terms  of  dipole  moments 
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Figure  16  shows  6  vortex  paths  near  a  Joukowsky  profile.  The  parameters 
were  chosen  in  such  a  way  that  the  profile  approximates  a  NACA  0012  as  has 
been  used  in  the  experiments.  The  coefficient  of  circulation  is  1. 


Obermeier  used  matched  asymptotic  expansions  to  fit  an  acoustic  field  to 
the  hydrodynamic  far  field.  A  detailed  description  of  this  method  can  be 
found  in  [1,2,3,10,11]. 


Using  nondimensional  coordinates  for  the  inner  (index  i)  and  outer  (in¬ 
dex  a)  range 
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and  transforming  to  the  coordinates 
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one  gets  the  following  result  for  the  sound  pressure: 
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being  the  Hankel  function  of  first  type  and  order  one,  which  is  the 
solution  for  a  dipol-like  cylindrical  sound  field.  Besides  a  constant  fac¬ 
tor,  F 2 )  is  the  Fourier  transform  of  4>  (X  ) 


3.  Experimental  Facilities 


3.1.  Duct,  Vortex  Generators  and  Profiles 


The  transonic  duct  (Fig.  2)  is  a  vacuum  type  wind  tunnel  with  two  vacuum 
tanks  of  132  and  24  m  ^  volume  (12).  The  maximum  measuring  time  for  the 
profile  investigation  arrangement  is  10  s  at  Mach  1  and  25  s  at  Mach  0.3.  A 
new  experiment  can  be  performed  every  4  minutes  (this  is  the  time  the  pump 
needs  to  evacuate  the  tanks). 

Air  is  sucked  in  from  the  laboratory  or  an  air  drying  unit.  The  velocity  in 
the  measuring  chamber  is  controlled  by  an  adjustable  diffuser.  In  the  smal¬ 
lest  cross  section  of  the  diffuser  the  air  reaches  Mach  1 ,  so  there  is  con¬ 
stant  velocity  throughout  the  measuring  time.  The  velocity  of  the  air  can 
be  chosen  from  v=0  to  v=600  m/s  if  a  Laval  nozzle  is  used. 

The  test  section  of  the  duct  (Fig.  6)  is  800  mm  long,  330  mm  high  and 
100  mm  wide.  The  upper  and  the  lower  side  of  the  chamber  have  slots  to  pro¬ 
vide  good  acoustical  damping.  The  vortex  generators  are  installed  in  the 
200  mm  front  section  (6a,c,e)  or  at  locations  nearer  to  the  window  (6b, d). 
The  profile  can  be  mounted  between  the  windows  downstream  of  the  genera¬ 
tors.  The  windows  have  a  diameter  of  230  mm  and  have  interferometric  qual¬ 
ity.  Instead  of  the  windows  a  pressure  measuring  plate’  with  an  array  of 
holes  at  15  mm  distance  can  be  used.  The  signals  of  up  to  16  Kulite  trans¬ 
ducers  can  be  digitized  simultaneously.  Thus  accurate  pressure  field  re¬ 
cordings  can  be  made. 

TV/o  profiles  are  used  for  the  investigations.  They  have  a  NACA  0012  shape 
with  120  mm  and  60  mm  chord  length,  respectively.  For  the  investigations 
they  are  placed  between  the  windows  at  different  angles  of  attack. 

Different  kinds  of  bluff  bodies  are  used  as  vortex  generators.  Their  wake 
is  a  Karman  vortex  street  consisting  of  strong  vortices  of  opposite  sense 
of  direction.  The  frequency  of  the  vortices  depends  on  the  size  and  shape 
of  the  bluff  body  (a  circular  cylinder  has  a  nearly  constant  Strouhal 
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number  of  0.2;  i.e.  the  distance  between  vortices  of  equal  sense  of  direc¬ 
tion  is  5  diameters).  Several  bluff  body  shapes  were  investigated,  it 
turned  out  that  a  square  cylinder  has  the  highest  amplitudes  and  best  con¬ 
stancy  of  Strouhal  frequencies  over  the  Reynolds  number  range  from  0.1  up 
to  0.5  millions. 

3.2.  Shock  Tube 


Like  other  authors  [13]  we  have  used  a  shock  tube  to  generate  starting  vor¬ 
tices.  The  shock  tube  consists  of  a  high  pressure  section  of  2100  mm  length 
and  a  low  pressure  section  of  3700  mm  (fig.  7).  The  width  of  the  tube  is 
100  mm  and  the  height  is  330  mm  (fig.  8).  The  distance  from  the  membrane  to 
the  measuring  chamber  divided  by  width  and  height  of  the  tube  give  length 
to  width  and  length  to  height  ratios  of  27  and  8. 

The  bottom  and  the  top  of  the  low  pressure  section  are  rounded  with  a  dia¬ 
meter  of  100  mm  to  prevent  shock  disturbances  (fig.  8).  The  measuring  cham¬ 
ber  is  the  same  as  in  the  wind  tunnel  experiment  (fig.  9).  The  same  measur¬ 
ing  devices  can  be  used.  ljae  low  pressure  section  is  evacuated  by  the  wind 
tunnel  vacuum  tank  to  a  vacuum  up  to  100  mbar.  The  high  pressure  section  is 
pressurized  up  to  1500  mbar.  Pressure  ratios  from  1  up  to  15  can  be  ob¬ 
tained  . 

Figure  8  shows  the  measuring  chamber.  A  GOE  508  profile  of  120  mm  chord- 
length  is  mounted  before  the  window  area.  The  GOE  508  profile  is  a  high 
lift  profile  with  flat  underside.  It  is  mounted  at  an  angle  of  7.5  degrees. 
Its  stationary  lift  coefficient  is  1.25.  5  mm  of  the  trailing  edge  of  the 
profile  are  visible  in  the  viewing  field.  The  60  mm  NACA  0012  profile  is 
mounted  38  mm  apart  at  the  same  height  as  the  trailing  edge  of  the  asym¬ 
metric  profile. 

3.3.  Computer  Installation 


The  computer  is  a  PDP  11/34.  It  has  connections  to  the  image  processing 
computer,  also  a  PDP  11/34,  to  the  institute  computer,  a  VAX,  and  to  the 
central  computer  of  Gottingen  University  and  Max-Planck-Institute,  a  Uni- 
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vac  1108.  Figure  4  shows  the  configuration  of  the  PDP  and  the  periphery 
[14). 

On  the  right  are  the  data  aquisition  components.  A  Camac  crate  controller 
(CA1 1 -FP , Digital  Equipment)  is  the  interface  to  16  transient  recorders  (Le 
Croy),  a  real  time  clock  and  a  digital  I/O  port.  The  transient  recorders 
have  8  K  storage  capacity  of  10  bit  data  at  data  aquisition  rates  of  up  to 
1  MHz.  They  are  used  to  record  the  pressure  data. 

An  IEC-bus  interface  is  used  to  record  the  static  pressure  with  the  help  of 
a  digital  voltmeter. 

The  third  I/O  component  is  the  controller  (LPA11-K,  Digital  Equipment)  of  a 
12  bit  analog/digital  converter.  Data  rates  up  to  50  KHz  can  be  stored  di¬ 
rectly  on  the  disc. 

The  storage  of  measured  data  is  performed  with  a  160  Megabyte  disc  (System 
Industries).  Long  term  storage  is  done  with  a  magtape. 

Fast  data  processing  is  done  with  the  array  processor  (Floating  Point  Sys¬ 
tems)  .  It  is  well  suited  for  fast  Fourier  transforms,  a  FFT  of  1024  points 
is  done  in  3  ms. 

For  visualisation  of  data  there  graphic  terminals  (VT100  and  Tektronix)  and 
a  hard  copy  unit,  connected  to  the  computer  by  a  terminal  multiplexer. 

The  interface  to  the  fast  data  line  to  the  image  processing  computer  is  a 
DRU-11C  (Digital  Equipment).  A  DMC-1 1  is  the  interface  to  the  data  line  to 
the  VAX  and  Univac.  The  VAX  is  used  for  program  development  and  other  com¬ 
puting  purposes.  The  Univac  is  needed  for  DISSPLA  plots. 

3.4.  Optical  Facility 

A  Zeiss  Mach-Zehnder  interferometer  can  be  used  for  optical  investigations 
(fig.  3).  The  channel  width  is  100  mm.  Vacuum  in  the  measuring  chamber  cor¬ 
responds  to  48  whole  fringes  in  interferograms.  The  light  source  is  a  high 
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frequency  spark  light  source.  Movies 
ra  at  10  Khz  or  with  a  Fastax  camera 


can  be  made  with  a  rotating  drum  came- 
at  speeds  up  to  8  Khz. 


4.  Testing  of  the  Experimental  Arrangement 
4.1.  Vortex  Street  Generators 


Vortex  streets  should  be  two-dimensional,  should  have  a  constant  Strouhal 
number  and  must  consist  of  strong  potential  vortices  of  the  right  spacing 
to  be  of  use  in  our  experiment.  The  interaction  of  vortices  among  them¬ 
selves  should  be  less  than  the  interaction  of  profile  and  vortex.  This 
means  that  the  quotient  length/(2*Strouhal  number),  the  spacing  of  the  vor¬ 
tices,  should  be  at  least  of  the  same  magnitude  as  the  profile  length.  The 
40  mm  cylinder  has  a  Strouhal  number  of  0.2  and  shows  a  vortex  spacing  of 
40/(2*0.21=100  mm.  Although  the  aspect  ratio  of  this  cylinder  is  very 
small,  the  vortex  street  is  quite  stable  and  has  only  a  single  main  fre¬ 
quency.  Figure  10a  shows  the  frequency  at  different  pressure  ratios.  The 
Strouhal  number  is  nearly  constant  but  there  is  some  influence  of  resonance 
of  the  duct  geometry. 

The  40  mm  cylinder  and  the  120  mm  profile  were  taken  for  the  first  vortex- 
profile  interaction  tests.  Figure  10b  shows  the  arrangement  working  with  an 
even  higher  amplitude  than  the  pure  cylinder. 

For  tests  with  the  60  mm  profile  several  bluff  body  shapes  were  investi¬ 
gated  (fig.  5).  It  turned  out  that  a  square  cylinder  has  the  highest  ampli¬ 
tudes  and  best  consistency  of  Strouhal  frequencies  over  the  Reynolds  number 
range.  The  20  mm  cylinder  does  not  work  as  well  as  the  bigger  one,  perhaps 
because  of  a  different  aspect  ratio. 
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4.2.  Shock  Tube 


The  shock  tube  setup  uses  the  same  measurement  equipment  as  the  duct.  The 
measuring  chamber  is  the  same.  TTie  length  of  the  tube  is  5800  mm.  It  is 
sufficient  to  provide  measuring  times  of  5  ms  minimum  for  the  Mach 
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number  range  from  0.2  to  0.8.  The  undisturbed  trajectory  length  of  the 
vortex  is  at  least  400  mm. 

Because  of  the  height  of  the  chamber  there  are  restrictions  in  the  attain¬ 
able  length  to  height  ratio  of  the  tube.  The  length  from  the  membrane  to 
the  windows  is  2600  mm.  Thus  the  length/height  ratio  is  only  8.  As  mterfe- 
rograms  prove,  this  small  ratio  is  enough  to  develop  the  shock.  A  problem 
is  the  rectangle  geometry  of  the  tube.  The  bursting  of  the  rectangular  mem¬ 
brane  generates  a  circular  wave.  The  reflections  from  the  walls  are  oblique 
waves.  These  waves  steepen  up  and  form  oblique  shocks.  To  prevent  these 
disturbances  it  was  tried  to  change  the  bursting  process  of  the  membrane  by 
taking  other  materials  and  destroying  the  membrane  by  high  voltage.  But 
this  did  not  work  very  well  and  it  was  found  that  round  shapes  of  the  top 
and  the  bottom  of  the  tube  prevent  uniform  reflections  of  the  oblique 
disturbances . 


5.  Computer-based  interferogram  evaluation 
5.1.  Digitizing  and  preprocessing 
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Interferograms  are  recorded  photographically  by  taking  single  frames  or 
taking  series  of  interferograms  on  16  mm  film  using  high  speed  pickup  tech¬ 
niques.  Digitizing  of  these  interferograms  is  done  by  an  image  processing 
system  (IP  512,  Imaging  Technology),  equipped  with  2  frame  buffers  of  512x 
512x8  bits  capacity  each.  The  image  processor  is  connected  to  a  mini -compu¬ 
ter.  The  system  further  is  equipped  with  a  digitizing  tablet  and  a  graphic 
terminal  to  enable  manual  input  of  graphic  data.  Also,  a  computer-control¬ 
led  film  projector  is  used  when  consecutive  frames  of  a  series  of  interfe¬ 
rograms  are  to  be  digitized.  A  detailed  description  of  the  system  is  given 
in  (15,16  1.  The  digitized  TV-lines  are  fed  to  the  computer  and  are  proces¬ 
sed  sequentially  or  stored  on  disk  for  subsequent  processing.  Depending  on 
the  quality  of  the  interferograms,  different  enhancement  procedures  can  be 
used  to  increase  the  signal-to-noise  ratio,  to  correct  for  uneven  illumi¬ 
nation  or  to  increase  the  contrast  of  the  fringes.  These  preprocessing 
routines  can  be  done  on  the  image  processor,  too. 
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5.2.  Fringe  segmentation 

Fringe  segmentation  is  performed  in  two  steps.  First  the  gray  level  fringes 
are  converted  to  binary  levels  by  using  a  fixed  or  floating  threshold.  In 
the  next  step  the  left-  und  right-edge  points  of  the  fringes  are  collected 
and  stored  in  a  polygonal  data  structure  using  a  sequential  tracking  algor¬ 
ithm.  To  reduce  the  amount  of  data  to  store,  a  redundancy  reduction  process 
approximates  the  actual  polygon  by  a  subset  of  the  vertices  of  the  original 
polygon  within  a  given  range  of  tolerance. 

If  fringes  leave  the  field  of  view,  or  if  a  background  object  resides  in¬ 
side  i_ne  fringe  field  (e.g.,  an  airfoil),  the  visible  area  of  the  fringes 
can  be  handled  by  using  the  points  on  the  boundary  as  edge  points  of  the 
fringes.  The  corresponding  polygons  of  course  are  not  to  be  along  the  boun¬ 
dary  of  the  background  object,  since  they  may  have  different  fringe  orders. 
In  order  to  establish  the  boundary  test  in  a  quick,  easy  and  robust  manner, 
it  is  not  desirable  to  derive  the  boundary  information  from  the  fringe- 
field  itself.  Instead,  the  geometry  of  the  test  section  is  used  to  generate 
a  binary-valued  mask  which  is  compared  with  the  actual  interferogram  pixel 
by  pixel  while  performing  the  fringe  extraction  process. 

Fringe  spacing  may  be  very  small  at  some  locations  in  the  fringe  fields,  as 
is  the  case  for  instance  inside  the  boundary  layer,  inside  a  shock  wave  or 
near  separation  lines.  If  the  resolution  of  the  digital  system  is  exceeded 
at  these  locations,  disconnected  or  falsely  connected  fringes  may  occur. 
This  results  in  the  fact  that  some  of  the  polygon  segments,  representing 
different  fringe  orders,  may  be  linked  together.  In  order  to  improve  the 
numbering  process,  most  of  the  disconnections  can  be  removed  by  an  analysis 
of  some  geometrical  parameters  of  the  polygons.  These  geometrical  parame¬ 
ters  are  the  shape  feature  (the  circularity  defined  by  the  enclosed  area 
divided  by  the  perimeter  squared),  the  ratio  of  the  distance  of  the  polygon 
endpoints  to  the  length  of  the  polygon  line,  and  the  angles  between  the  po¬ 
lygon  segments. 

Those  lines  which  are  suspected  to  have  disconnections  are  cut  by  this  pro¬ 
cess  at  locations  where  the  polygon  is  folded  or  where  it  has  sharp  edges. 


The  remaining  polygons  having  disconnections  not  detected  by  this  cutting 
process  have  to  be  handled  by  the  numbering  process. 

Some  consideration  and  programming  work  has  been  done  on  the  problem  of 
getting  numbers  for  the  fringe  lines  of  an  interference  fringe  field  by 
using  the  known  numbers  of  a  similar  fringe  field.  The  matter  is  somewhat 
more  complicated,  as  some  of  the  lines  of  the  actual  fringe  field  may  have 
disconnections,  and  therefore  it  may  be  impossible  to  get  a  unique  fringe 
number  for  those  lines. 

Tbe  idea  of  a  numbering  scheme  applicable  to  series  of  interferograms  is  to 
get  additional  information  from  the  fringe  field  numbered  previously.  This 
may  only  be  possible  when  the  fringe  locations  vary  only  by  a  small  frac¬ 
tion  of  the  fringe  spacing  at  almost  any  locations  of  the  test  section.  If 
the  interframe  time  between  consecutive  interferograms  is  choosen  appro¬ 
priately,  this  condition  is  satisfied  in  the  case  of  the  profile  flow 
investigation.  The  numbering  of  a  series  of  interferograms  starts  with  the 
setting  of  the  numbers  of  the  first  fringe  field  by  hand.  The  fringe  lines 
are  superimposed  by  a  set  of  so  called  test-lines  at  which  a  set  of  order 
number  functions  is  defined  by  a  rational  spline  approximation.  This  set  of 
spline  functions  is  fitted  to  the  fringe  lines  of  the  following  interfero- 
gram.  The  numbers  of  those  fringes  overlapping  uniquely  at  the  test-lines 
are  accepted  directly,  the  other  numbers  are  derived  from  the  spline  func¬ 
tions.  An  older  version  of  the  program  uses  a  set  of  straight  lines  running 
in  the  x-  and  y-direction,  while  a  newer  version  uses  a  set  of  polygons  as 
test-lines.  The  advantage  of  the  polygonal  test-lines  is  that  they  are  more 
suitably  adaptable  to  the  global  course  of  the  fringe  lines. 

In  the  case  where  some  disconnected  lines  remain  after  preprocessing  a 
unique  number  cannot  be  assigned  to  these  lines  by  the  numbering  scheme.  To 
decide  where  these  lines  have  to  be  divided  into  two  or  more  segments  -  an 
older  version  of  the  program  uses  the  angles  between  subsequent  polygon 
segments  -  the  deviations  from  a  local  approximation  of  the  fringe  order 
function  is  used.  Another  advantage  of  comparing  fringe  fields  by  using 
local  polynomial  approximations  is  the  fact  that  fringe  polygons,  having  no 
intersections  with  the  test-lines,  could  be  processed  using  the  same  algo- 


nthm,  leading  to  a  more  stable  behaviour  of  this  numbering  process.  As  the 
position  of  the  airfoil  is  not  fixed  relative  to  the  video  coordinates'  of 
the  scanning  equipment  (due  to  imprecise  location  of  successive  frames  on 
the  film  inside  the  high  speed  camera  and  other  causes),  software  had  to  be 
developed  to  achieve  automatic  adjustment  of  image  coordinates  and  masking 
areas.  This  is  especially  important  when  large  series  of  interferograms  are 
to  be  processed  automatically. 

5.3.  Polynomial  approximation 

The  fringe  order  function  is  defined  at  a  set  of  contour  lines,  but  most  of 
the  mathematical  transformations  to  follow  require  the  interpolation  of 
fractional  fringe  order  numbers.  Some  methods  to  interpolate  between  ran¬ 
domly  scattered  points  on  a  surface  are  discussed  in  the  literature  (17). 
We  have  used  a  local  distance  weighted  polynomial  least  squares  approxima¬ 
tion,  best  suited  in  regard  to  computation  time  and  numerical  representa¬ 
tion.  This  process  computes  the  coefficients  of  two-dimensional  polynomials 
of  second  order  at  the  mesh  points  of  a  rectangular  grid.  To  calculate  the 
interference  order  inside  the  meshes,  the  four  polynomials  at  adjacent  cor¬ 
ners  of  the  mesh  are  evaluated  and  weighted  proportionally  to  their  dis¬ 
tance  to  the  point  of  interest. 

Another  method  is  to  compute  a  two-dimensional  spline  approximation,  using 
some  of  the  previously  computed  polynomial  coefficients  to  get  a  smooth 
surface  having  continuous  first  and  second  derivatives.  The  density  func¬ 
tion  may  easily  be  computed  by  a  linear  transformation  of  the  fringe  order 
values,  because  in  the  actual  experimental  set-up  the  fringes  are  lines  of 
constant  density  (infinite  fringe  case).  In  the  case  of  finite  fringe 
fields  a  subtraction  of  the  overlayed  fringe  pattern  could  be  achieved  by  a 
modification  of  some  of  the  polynomial  coefficients. 

Routines  have  been  developed  to  add  or  subtract  the  spline  approximations 
describing  different  fringe  fields  recorded  under  equivalent  physical  con¬ 
ditions.  Thus  many  fields  may  be  averaged  and  compared  with  theoretical 
calculations.  This  method  has  been  applied  to  areas  where  the  flow  field  of 
interest  had  been  disturbed  by  turbulence.  When  the  basic  flow  field  chan- 


ges  slowly  as  compared  with  typical  time  scales  of  turbulence,  and  when  the 
photographic  frequency  is  high  enough,  averaging  can  be  done  over  consecu¬ 
tive  frames  to  eliminate  turbulence.  In  principle,  turbulence  could  be  in¬ 
vestigated,  too,  by  comparing  single  fields  with  an  average  (figure  35). 

5.4.  Localization  of  the  vortex  core 


In  the  study  of  the  vortex-profile  interaction  process  a  point  of  interest 
is  the  knowledge  of  the  vortex  traces  under  different  positions  of  the  pro¬ 
file,  to  see  whether  a  theoretical  description  of  these  traces  is  in  accor¬ 
dance  with  experimental  results.  Some  program  development  has  been  done  to 
localize  the  typical  fringe  pattern  resulting  from  the  appearance  of  the 
vortex.  The  vortex  core  is  computed  using  the  centroid  of  the  enclosed  area 
of  the  innermost  of  these  fringe  lines. 

5.5.  Evaluation  of  interferograms 


As  an  example  of  digital  interferogram  analysis  two  processed  fringe  fields 
of  a  profile  flow  with  and  without  a  vortex  present  are  shown  in  Figures 
11a  and  12a.  The  flow  velocities  are  280/s  and  100  m/s,  respectively,  at  an 
angle  of  attack  of  0  degrees.  The  interferograms  were  digitized  with  a  re¬ 
solution  of  512  by  256  pixels.  The  fringe  extraction  process  described  in 
chapter  5.2  was  used  to  obtain  the  polygonal  fringe  fields  shown  m  Figures 
11b  and  12b.  The  boundaries  of  the  test  section  as  well  as  those  of  the 
airfoil  are  input  by  hand,  using  the  digitizing  tablet.  These  lines  are 
used  to  prepare  a  mask  file,  needed  by  the  fringe  extraction  program  to 
avoid  connections  of  fringe  lines  along  the  boundaries  of  background  ob¬ 
jects.  Due  to  the  restricted  resolution  of  the  digital  system  there  ap¬ 
peared  some  falsely  connected  lines  around  the  leading  edge  of  the  airfoil, 
where  the  fringe  spacing  is  very  small.  These  locations  could  be  detected 
and  the  falsely  connected  lines  could  be  cut  by  application  of  the  afore¬ 
mentioned  algorithm,  using  the  geometrical  parameters  of  the  lines. 

The  lines  are  numbered  corresponding  to  the  interference  order,  which  in 
this  case  is  done  by  hand,  because  the  automatic  numbering  algorithm, 
applicable  to  these  types  of  fringe  fields,  needs  a  numbered  field  to  start 
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the  numbering  process  of  consecutive  frames.  Hie  density  function  was 
derived  from  the  fringe  field  by  approximation  of  polynomials  of  second 
order  at  the  mesh  points  of  a  rectangular  grid  (shown  in  figs.  11b  and 
12b)  as  described  in  chapter  5.2.  Tto  achieve  a  smooth  surface,  the  inter¬ 
ference  order  function  was  calculated  from  these  polynomials  using  a  bicu¬ 
bic  spline  interpolation  procedure.  Two  density  profiles  in  x-direction  at 
y  =  0.2  are  plotted  in  Figs.  11c  and  12c.  Tbe  relative  density  as  a  func¬ 
tion  of  the  fringe  order  is  given  in  this  case  as 

1  pc  Tn  X  N  (  x ,  y ) 


where  o/o  is  the  relative  density,  n  (T,  p)  is  the  refractive  index  of  air 
at  temperature  T  and  pressure  p,  T  and  p  are  the  temperature  and  pressure 
at  rest,  is  the  wavelength  of  light,  h  the  depth  of  the  test  chamber  and 
N(x,y)  the  fringe  order  function.  Note  that  the  order  function  is  divided 
by  two,  because  each  fringe  is  represented  by  two  lines  (a  "left"  and  a 
"right"  one).  The  above  equation,  evaluateo  with  the  actual  dimensions 
valid  in  the  tests,  yields 


o/o  =  1  =  0.01058  *  N  (x,y) 

As  can  be  seen  from  Fig.  17,  the  lowest  density  appears  inside  the  vortex 
core  with  a  value  of  0.92  relative  to  the  density  at  the  stagnation  point, 
where  the  fringe  order  was  set  equal  to  zero. 


6.  Experimental  results  of  vortex-street-prof lie  interaction 


6.1.  Pressure  fields 

Measurements  were  performed  with  the  40  mm  cylinder  with  and  without  the 
120  mm  profile.  Tbe  pressure  ratio  range  was  from  0.999  to  0.6,  that  is 
from  13  m/s  to  290  m/s. 


\-\Al 
■  c  o 


As  expected,  the  pressure  field  is  antisymmetric  when  vortex  generator  and 
profile  are  aligned  in  the  cord  plane.  The  difference  from  antisymmetry'  is 
less  than  5  degrees  phase  angle  in  Fourier  spectra.  Therefore  it  is  suffi- 
|  cient  to  have  pressure  transducers  only  in  the  upper  half  of  the  measuring 

chamber . 

Figure  6a  shows  the  pressure  transducer  positions  used  in  the  measurements 
which  are  described  in  this  report.  Figure  10  is  done  with  pressure  trans- 
!  ducer  5.  Figure  13a  and  b  are  typical  examples  of  pressure  recordings.  At  a 

pressure  ratio  of  0.83  transducer  2  measures  an  underpressure  of  more  than 
;  200  mbar.  The  vortices  differ  in  strength  and  do  not  have  the  same  path,  so 

^  the  extreme  underpressure  peaks  only  happen  from  time  to  time.  There  is  a 

*  diffusion  of  vortex  strength  and  the  peaks  at  transducers  8,10  and  12  do 

;  not  have  such  high  amplitudes  any  more.  Signals  at  transducers  5,6,7  and  9 

■  (above  profile)  do  not  have  single  peaks  any  more  but  a  more  sinusoidal 

"  shape . 

r 

;  Depending  on  resonance  the  peaks  of  frequency  recordings  are  more  or  less 

;  sharp.  But  there  is  always  a  Strouhal  frequency  far  above  the  wall  turbu- 

* 

•  lence  level  at  nearly  all  measuring  points.  Figures  14  to  1 7  are  typical 

|  examples  of  spectra.  The  spectra  are  done  with  0.8  s  long  pressure  record¬ 

ings  at  10  Khz  sample  frequency. 


At  the  same  pressure  ratios  vortex  street  and  profile  vortex  interaction 
have  the  same  Strouhal  frequency.  The  vortex  generation  process  seems  not 
to  be  influenced  by  the  profile.  The  profile  measurements  show  only  a 
slightly  better  consistancy  of  Strouhal  frequency.  The  amplitudes  in  the 
vicinity  of  the  cylinder  are  of  the  same  order  of  magnitude.  The  amplitudes 
measured  above  the  profile  (transducers  3, 5, 6, 7, 9)  are  greater  in  the  pro¬ 
file  case.  At  any  pressure  ratio  the  phase  differences  between  the  measur¬ 
ing  points  are  entirely  different.  The  hydrodynamic  field  as  well  as  the 
sound  field  are  influenced  strongly  by  the  profile. 


The  first  pair  of  spectra  (Fig.  14)  is  done  at  the  high  pressure  ratio  of 
0.998.  The  velocity  is  about  19  m/s.  Sound  generation  is  inefficient  at  low 
velocities.  The  pressure  field  is  a  relatively  undisturbed  hydrodynamic 
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field.  The  Strouhal  frequency  is  of  the  same  magnitude,  but  the  vortex 
street  signals  are  more  noisy.  Very  obvious  is  the  difference  at  transducer 
1.  Only  transducer  8  of  the  profile  case  shows  a  noisy  character,  because 
it  is  influenced  by  the  profile  wake.  It  is  interesting  that  the  phase  dif¬ 
ferences  between  the  same  transducer  pairs  are  entirely  different  even  in 
the  pure  hydrodynamic  case. 

The  second  pair  (fig.  15)  has  a  sharp  dominant  frequency  with  the  profile 
and  several  frequencies  with  the  vortex  street,  due  to  beats  and  frequency 
jumps.  Figure  16  is  an  opposite  example.  The  vortex  street  shows  the  more 
narrow  frequency  band,  whereas  the  profile  has  got  a  large  number  of  diffe¬ 
rent  Strouhal  frequencies.  Surprisingly,  there  is  the  same  dominant  fre¬ 
quency  at  all  measuring  points.. 

Figure  1 7  is  a  resonance  case  with  extremly  high  amplitudes.  There  is  the 
same  frequency  at  all  measuring  points  and  nearly  the  same  frequency  in 
both  plots.  (The  first  harmonic  has  a  higher  amplitude  at  transducer  1  and 
8,  but  the  fundamental  is  the  same).  In  both  cases  the  highest  amplitudes 
are  at  transducers  2  and  3.  The  transducers  above  the  profile  (3, 5,6, 7,9) 
have  higher  amplitudes  with  the  profile.  The  greatest  difference  is  in  the 
first  harmonic  (fig.  18).  There  is  a  harmonic  at  transducers  1,2,3  in  both 
cases,  but  the  profile  shows  higher  amplitudes  and  higher  quality  factor. 
Entirely  different  are  the  transducer  signals  above  the  profile  (5, 6, 7, 9) 
where  the  vortex  street  shows  noise  and  the  profile  shows  the  first 
harmonic  with  high  amplitudes  and  high  quality  factor.  This  might  be  a  hint 
that  a  new  phenomenon  is  occurring  with  strong  vortices  at  high  Reynolds 
numbers. 

6.2.  Phenomena  of  interaction 

The  boundary  layer  plays  an  important  role  in  vortex  airfoil  interaction. 
The  closer  the  vortex  gets  to  the  airfoil,  the  stronger  is  the  effect  on 
the  boundary  layer.  Figure  19  is  an  overview  on  the  different  interaction 
phenomena.  (The  flow  comes  from  the  left  side,  the  sense  of  direction  of 
the  vortex  is  counterclockwise).  The  interaction  is  shown  in  4  time  steps. 
Depending  on  the  distance  to  the  airfoil  axis,  we  can  find  eight  (a-h) 
principally  different  cases.  The  phenomena  are  very  much  alike  for  Mach 


numbers  from  0.2  to  0.8  and  circulation  coefficients  from  0.3  to  1.5.  The 
case  f  happens  only  at  circulation  coefficients  above  0.2.  Case  c  needs  a 
circulation  coefficient  above  0.5  or  a  high  Mach  number. 

Cases  a  and  h:  There  is  nearly  no  interaction  with  the  airfoil  boundary 
layer  at  a  greater  distance  from  the  airfoil  axis,  that  is  more  than  one 
third  of  airfoil  length  at  a  circulation  coefficient  of  0.5.  Only  the  stag¬ 
nation  point  is  moving  downwards  and  a  tiny  suction  region  is  generated  on 
the  upper  side  of  the  airfoil.  After  the  vortex  has  passed  the  leading  edge 
there  are  no  strong  deviations  from  the  stationary  flow  around  the  airfoil. 

Cases  b  and  g:  When  the  vortex  comes  closer  to  the  airfoil  the  interfero- 
gram  shows  a  small  thickening  in  boundary  layer  between  vortex  and  airfoil 
boundary.  The  thickening  is  generated  at  about  one  quarter  of  airfoil  cord- 
length.  It  is  travelling  with  the  vortex  and  it  is  slowly  growing.  This 
thickening  can  be  the  source  of  a  secondary  vortex  when  the  interaction  is 
strong  enough.  This  secondary  vortex  is  generated  downstream  of  the  primary- 
vortex  if  the  vortex  is  above  the  airfoil  and  upstream  if  the  vortex  is  be¬ 
low  the  airfoi 1 . 

Cases  c  and  f:  If  the  vortex  is  stronger  there  is  a  separation  in  the  suc¬ 
tion  peak  which  becomes  a  secondary'  vortex.  Case  f,  where  the  vortex  passt-s 
under  the  airfoil  occurs  only  at  circulation  coefficients  above  2.0.  At 
very  high  circulation  coefficients  of  about  4.0  there  is  separation  at  the 
suction  peak  if  the  vortex  is  still  more  than  half  an  airfoil  length  away 
from  the  leading  edge. 

Cases  d  and  e:  If  the  vortex  hits  the  leading  edge  directly,  it  vanishes 
totally  (Case  e).  If  the  vortex  passes  very  close  to  the  profile  the  vortex 
looses  a  lot  of  its  circulation  and  only  a  very  small  turbulent  spot  is 
lef t  (Case  d ) . 

The  figures  20,21,22  and  23  are  examples  for  the  cases  b,  c,  e  and  f.  'Hie 
Figures  20,21  and  22  are  related  to  the  shock  tube  experiment  while  fig.  12 
comes  from  the  stationary  wind  tunnel  experiment.  Figure  20  contains  ten 
frames.  A  vortex  is  coming  from  the  left  side  and,  while  passing  the  an- 
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foil,  it  generates  a  secondary  vortex.  Figure  21  and  22  are  experiments  at 
the  high  Mach  number  of  0.8.  In  the  first  of  the  5  frames  the  shock  wave 
and  the  vortex  are  at  a  close  distance.  Before  the  vortex  gets  into  contact 
with  the  airfoil  it  passes  the  leading  edge  reflection  of  the  shock.  When 
the  vortex  gets  close  to  the  leading  edge  there  is  a  strong  interaction.  In 
both  cases  the  vortex  nearly  vanishes.  But  this  strong  interaction  gener¬ 
ates  a  strong  acoustic  wave.  Figure  23  shows  an  interaction  of  the  airfoil 
with  a  vortex  of  the  high  circulation  coefficient  of  4.  There  appears  a 
strong  secondary  vortex  generated  in  the  suction  region  of  the  airfoil. 

6.3.  Experimentally  determined  vortex  paths 


According  to  theoretical  calculations  as  described  above  and  by  other  au¬ 
thors  [2,11,6)  the  generation  of  sound  depends  on  the  path  of  the  vortex. 
Therefore,  precise  knowledge  of  these  paths  is  necessary  for  testing  these 
theories.  Figures  24,25,26  show  paths  of  vortices  under  different  condi¬ 
tions  in  comparison  with  streamlines  of  incompressible  flow  around  a  Jou- 
kowsky  profile.  Short  streamlines  near  the  trailing  edge  indicate  secondary 
vortices . 

Far  away  upstream  the  vortices  move  uniformly  with  the  flow.  Near  the  lead¬ 
ing  edge  they  deviate  from  the  streamlines  in  an  upward  direction.  This  de¬ 
viation  is  larger  at  small  Mach  numbers  then  at  large  pnes.  This  agrees 
with  a  prediction  by  Lindblad  (181. 

The  behaviour  of  the  vortex  after  passing  the  leading  edge  depends  on  the 
distance  from  the  profile.  If  the  vortex  hits  the  leading  edge  precisely, 
it  disappears.  When  the  distance  is  small,  a  strong  interaction  with  the 
boundary  layer  takes  place  and  secondary  vortices  are  generated.  This  rein¬ 
forces  the  upward  movement  of  the  vortex  above  the  upper  face  of  the  pro¬ 
file;  this  movement  continues  after  the  vortex  has  passed  the  trailing 
edge.  (An  exception  to  this  is  the  case  M  =  0.4,  Z  =  2.0,  where  the  secon¬ 
dary  vortex  is  stronger  than  the  residue  of  the  primary  one,  which  does 
rotate  around  the  former). 


The  appearance  of  a  secondary  vortex  enlarges  the  distance  between  vortex 
and  profile  under  the  lower  face,  too. 
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With  increasing  distance  from  the  profile  secondary  vortices  are  not  gener¬ 
ated  any  more.  The  vortex  reaches  its  uppermost  point  approximately  above 
the  centre  of  the  profile.  Then  it  moves  downwards  as  compared  with  the 
streamlines.  The  path  of  the  vortex  is  approximately  symmetrical  with  re¬ 
spect  to  the  mid-perpendicular  of  the  profile. 

The  vortex  paths  calculated  with  the  model  described  above  are  in  qualita¬ 
tive  agreement  with  the  theoretical  results  when  the  vortex-airfoil  sepa¬ 
ration  is  large.  However,  a  quantitative  comparison  shows  grave  disagree¬ 
ment.  The  first  and  third  plot  of  figure  27  show  a  close  contiguity  of 
theoretical  and  experimental  paths,  but  only  because  both  don't  deviate 
from  the  streamlines  very  much. 

Figure  28  shows  an  example  of  large  deviation  from  the  streamlines.  As  is 
the  case  in  the  other  examples,  the  theoretical  values  are  nearer  to  the 
streamlines  than  the  experimental  ones.  The  second  plot  in  Figure  27  shows 
a  comparison  between  experimental  data  and  theoretical  results  calculated 
without  taking  into  account  a  central  vortex  of  circulation.  (In  this  case 
the  circulation  around  the  profile  doesn’t  equal  zero  but  the  negative 
value  of  the  circulation  of  the  vortex.) 

Calculated  deviations  of  the  vortex  path  from  the  streamlines  are  larger 
than  in  the  model  with  a  central  vortex  and  larger  than  the  experimental 
values,  too.  Reality  seems  to  lie  somewhere  between  these  two  possibili¬ 
ties.  Circulation  around  the  airfoil  apparently  changes  when  a  vortex  is 
passing.  Therefore,  no  reasonable  agreement  can  be  achieved  when  the  model 
described  above  is  used. 

It  can  be  concluded  that  simple  models  of  potential  theory  describing  the 
passage  of  a  vortex  near  an  airfoil  will  meet  basic  difficulties.  To 
achieve  a  better  understanding,  precise  knowledge  of  the  influence  of  the 
boundary  layer  will  be  necessary.  A  further  difficulty  is  that  the  mecha¬ 
nism  of  varying  circulation  of  the  airfoil  is  not  yet  understood. 

6.4.  Generation  of  sound 
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Theory  predicts  the  generation  of  sound  with  wavelength  of  the  order  of 


magnitude  velocity  of  sound  (time  of  passage)  when  the  height  of  the  vortex 
above  the  horizontal  axis  of  the  airfoil  exceeds  0.1.  Tbe  wavelength  then 
is  larger  than  the  profile,  which  is  an  assumption  of  this  theory  [1,2).  As 
the  paths  of  the  vortices  are  not  predicted  accurately,  one  cannot  hope  for 
quantitatively  correct  results.  However,  some  qualitative  agreement  can  be 
expected. 

Waves  of  such  a  low  frequency  cannot  be  detected  in  interferograms,  but 
earlier  meassurements  indicated  the  appearance  of  low  frequency  sound  [19). 

In  addition  to  low  frequency  sound,  sound  waves  of  very  short  wavelength 
are  generated.  Shocks  appear  at  higher  Mach  numbers  (0.6,  0.8;  0.4  in  some 
special  cases)  these  shocks  emanate  from  the  leading  edge  of  the  profile. 
As  the  flow  near  the  leading  edge  does  not  depend  strongly  on  the  flow  at 
other  parts  of  the  airfoil,  the  wavelength  of  this  sound  presumably  depends 
on  the  radius  of  curvature  of  the  leading  edge  and  the  core  radius  of  the 
vortex . 

Figure  29  shows  4  series  of  interferograms  recorded  under  equal  conditions. 
The  vortex  directly  hits  the  leading  edge.  The  first  picture  of  every 
series  shows  the  shock  wave  and  its  reflection  off  the  leading  edge.  The 
interaction  of  the  vortex  with  the  reflected  cylindnal  wave  is  apparent  in 
the  second  picture  of  the  first  series.  As  a  consequence  of  this  interac¬ 
tion,  a  secondary  wave  behind  the  reflected  shock  appears  in  the  consecu¬ 
tive  pictures  which  is  inclined  to  the  right  by  approximately  10  -  15 
degrees  as  compared  with  the  reflected  shock.  The  second  pictures  of  the 
three  other  series  show  the  vortex  core  immediately  before  hitting  the 
leading  edge.  The  vortex  disappears  almost  completely;  only  weak  residues 
can  be  discerned  (3.  picture  of  2.  and  3.  series).  Instead,  a  strong  cylin¬ 
drical  wave  has  appeared,  the  centre  of  which  is  located  near  the  leading 
edge.  This  wave  steepens  to  form  a  shock  wave  (4.  and  5.  picture). 

A  binary  structure  is  characteristic  for  the  generated  wave  (3.  picture, 
4.  series).  The  length  of  the  wave  is  comparable  to  the  width  of  the  high 
velocity  region  immediately  behind  the  leading  edge  (2.  pictures).  This 
region  vanishes  too  fast  for  the  process  to  be  resolved  with  a  photographic 
frequency  of  10  kHz. 
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The  stagnation  point  jumps  to  the  front  which  results  in  a  pressure  jump  at 
the  leading  edge.  The  increase  of  pressure  at  the  upper  face  is  much  stron¬ 
ger  than  the  decrease  at  the  lower  one,  so  a  monopole-like  emission  is 
[  likely.  The  region  of  increasing  pressure  extends  from  the  leading  edge  to 

approximately  1/4  of  the  chordlength.  The  thickness  of  this  region  -  and 
therefore  a  characteristic  length  of  the  wavelength  of  generated  sound  -  is 
smaller  than  its  extension  in  the  direction  of  flow.  Therefore,  emission 
must  take  place  in  the  normal  direction.  Figure  30  shows  4  consecutive 
stages  of  propagation  of  the  generated  wave. 

If  one  approximates  the  wave  by  a  cylindrical  one,  the  centre  of  the  cylin¬ 
der  is  located  near  the  lower  face  and  moves  to  the  right  with  the  velocity 
of  the  flow.  An  angular  characteristic  has  been  plotted  for  the  circle  with 
its  centre  indicated  by  a  cross.  The  outer  circle  corresponds  to  a  density 
difference  of  one  fringe,  approximately  are  third  of  the  stagnation  densi¬ 
ty.  The  maximum  of  the  lobe  is  at  appr.  45  degrees  in  the  upper  left.  The 
maximal  amplitude  is  at  appr.  0.7  fringes,  corresponding  to  0.2  times  the 
stagnation  pressure.  As  the  infinite  fringe  mode  of  measurement  has  been 
used,  only  a  few  fringes  appear.  Therefore,  experimental  uncertainties  are 
quite  large.  The  maximum  amplitudes  of  the  four  waves  differ  by  about 
50  %.  The  form  of  interaction  depends  sensitively  on  the  height  of  the  vor¬ 
tex  above  the  horizontal  axis  of  the  airfoil  in  the  case  of  direct  hits. 

The  location  of  the  centre  of  the  wave  under  the  airfoil  indicates  that  the 
wave  is  emitted  off  an  extended  region,  so  the  centre  appears  to  be  situ¬ 
ated  behind  the  actual  region  of  generation. 

When  the  vortex  does  not  hit  the  leading  edge  directly,  more  complicated 
processes  take  place.  The  vortex  does  not  disappear  but  interacts  with  the 
low  pressure  region  of  the  airfoil.  Pressure  increases  at  the  upper  face  in 
this  case,  too,  but  the  region  of  maximal  pressure  change  cannot  be  located 
as  precisely  as  in  the  case  of  a  direct  hit. 


TVo  cases  of  a  vortex  passing  above  the  profile  are  shown  in  Fig.  31 .  The 
vortex  passes  far  away  from  the  airfoil  axis  in  the  series  at  the  left.  An 
interaction  of  type  c)  takes  place.  A  low  pressure  region  to  the  right 


of  the  vortex  can  be  discerned  on  picture  2.  A  secondary  vortex  and  a  com¬ 
pression  wave  become  apparent  (picture  3).  Tbe  other  series  is  more  similar 
to  the  first  case.  Although  the  vortex  does  not  disappear  completely,  the 
disappearance  of  the  low  pressure  region  is  analogous  to  the  case  of  a  di¬ 
rect  hit.  Tbe  wedge  of  low  pressure  downstream  above  the  vortex  which  ap¬ 
pears  in  conjunction  with  vortex  generation  can  be  seen  in  Figure  31  and 
should  not  be  confused  with  a  sound  wave. 

Figs.  32a  and  32b  show  a  cross  section  of  the  pressure  and  an  arxgular  char¬ 
acteristic.  The  waves  are  circular  in  both  cases,  their  center  lying  a  bit 
beneath  the  centre  of  the  airfoil. 

Ibe  amplitude  is  smaller  than  that  of  the  wave  in  32b  by  a  factor  of  about 
2,  their  radii  being  about  equal.  The  angular  characteristic  is  similar  in 
both  cases.  Because  of  the  stronger  circulation  of  the  vortex,  the  ampli¬ 
tude  is  larger  in  case  32b  than  in  that  of  30.  As  the  radius  of  the  wave  is 
larger  than  in  30,  too,  the  relation  between  circulation  and  amplitude  ap¬ 
pears  to  be  approximately  quadratic. 

Figs.  21,22  show  a  comparison  between  a  vortex  passing  over  and  under  the 
airfoil,  respectively.  Angular  characteristics  are  shown  in  Fig.  33. 

The  density  on  the  upper  face  of  the  airfoil  has  been  plotted  for  7  points 
in  Fig.  34.  Pressure  decreases  when  the  vortex  is  approaching  (example  at 
the  left);  there  is  a  steep  increase  of  pressure  after  the  passage  of  the 
vortex.  As  has  been  explained  above,  these  pressure  changes  cause  the  gen¬ 
eration  of  sound.  However,  the  temporal  resolution  of  the  measuring  equip¬ 
ment  is  not  yet  sufficient  to  resolve  the  process  of  generation  in  suffi- 
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Figure  8  ;  Front  v/ieu  of  the  shock  tube.  The  bottom  and 
the  top  of  the  tube  are  round  to  prevent  shock 
disturbances. 
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Figure  10:  Fourier  spectra  at  different  pressure  ratios. 
10a:  Vortex  street 

Transducer  5  Maximal  Amplitude:  14.17  [mbar] 
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10  b:  Profile  -  Vortex  Street  Interaction 

Transducer  5  Maximal  Amplitude:  49.92  [mbar] 
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Figure  13  a  and  13  b 
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40MM  CYLINDER  VORTE ^STREET  120MM  PROntE  IN  40MM  CYLVORTEXSTREET.12TR. 
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40MM  CYLINDER  VORTEX^TREET  120MM  PROFILE  IN  40MM  CYL.VORTEXiTREET.12TR. 
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Figure  18a  and  18b 


figure  20  Mach  0.2,  Z  =  0.5,  Re  =  290  000 
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Calculated  and  measured  paths. 

Upper  figure:  Mach  0.2, 

Re  =  290  000,  Z  =  0.5 

Middle  figure:  as  in  the  first  case  but  calculation  done 
without  the  term  corresponding  to  the  central  vortex. 
Lower  figure:  Mach  0.4,  Re  =  370  OOO,  Z  =  0.65 
Measured  curves  are  solid. 
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Figure  31 :  Generation  of  sound  by  vortices  passing  in  different  heights 

Mach  0.8,  Re  =  280  000,  Z  =  0.7. 

Left  series:  h  =  0.5;  right  series:  h  =  0.05. 


Figure  32a:  Density  along  a  line  through  the  center  of  the  wave;  angular 
diagram  for  31  (left  series).  Outer  circle  corresponds  to  one 
fringe  (one  third  of  stagnation  pressure). 
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files.  Thin  curves  show  the  vortex  paths.  Outer  circle 
responds  to  one  fringe  (one  third  of  stagnation  presstr 
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